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I. INTRODUCTION 



The Kondo problem describing the effect due to the exchange interaction between the magnetic impurity and the 
conduction electrons plays a very important role in condensed matter physics Q . Wilson g] developed a very powerful 
numerical renormalization group approach, and the model was also solved by the coordinate Bethe ansatz method 
which gives the specific heat and magnetization. More recently, a conformal field theory approach was developed 
by Affleck and Ludwig Q based on a work by Nozieres Q. In the conventional Kondo problem, the interaction 
between conduction electrons is discarded, due to the fact that the interacting electron system can be described as 
a Fermi liquid. Recently there has been substantial research devoted to the investigation of the theory of impurities 
coupled to Luttingcr liquids. Such a problem was first considered by Lee and Toner Q. By using the perturbative 
renormalization group theory they found that the Kondo temperature crosses from a generic power law dependence on 
the Kondo coupling constant to an exponential one in the infinite limit. Afterwards, a "poor man's" scaling procedure 
was carried out by Furusaki and Nagaosa |8| , who found a stable strong coupling fixed point for both antiferromagnetic 
and ferromagnetic cases. On the other hand, boundary conformal field theory predicts two types of critical behaviours, 
i.e., either a local Fermi liquid with standard low-temperature thermodynamics or the non- Fermi liquid observed by 
Furusaki and Nagaosa ||. However, in order to get a full picture about the critical behaviour of Kondo impurities 
coupled to Luttingcr liquids, some simple integrable models which allow exact solutions are desirable. 

Several integrable magnetic or nonmagnetic impurity problems describing impurities embedded in systems of corre- 
lated electrons have so far appeared in the literature. Among them are versions of the supersymmetric t — J model with 
impurities p- |l3| . Such an idea to incorporate an impurity into a closed chain dates back to Andrei and Johanesson 
(see also |lf||l6|]). However, the model thus constructed suffers from the lack of backward scattering and results 
in a very complicated Hamiltonian which is difficult to be justified on physical grounds. Therefore, as observed by 
Kane and Fisher pJ^ , it is advantageous to adopt open boundary conditions with the impurities situated at the ends 
of the chain when studying Kondo impurities coupled to integrable strongly correlated electron systems ]l8|~|20|] . 

In this paper, integrable Kondo impurities with arbitrary spin coupled to the one-dimensional t — J open chain 
are constructed following a formalism recently advocated in p(| . Our new input is to search for integrable boundary 
K matrices depending on the local magnetic moments of impurities, which arise as a nontrivial realization of the 
Z2-graded reflection equation (RE) algebras in a finite dimensional quantum space, which may be interpretated as 
an impurity Hilbert space. It should be emphasized that our new non-c-number boundary K matrices are highly 
nontrivial, in the sense that they can not be factorized into the product of a c-number boundary K matrix and the 
corresponding local monodromy matrices. The models we present are solved by means of the algebraic Bethe ansatz 
method and the Bethe ansatz equations are derived. 

The layout of this paper is the following. We begin by reviewing the Z2-gradcd boundary Quantum Inverse 
Scattering Method (QISM) as formulated in p5J|. We then introduce two integrable cases of the one-dimensional t — J 
model with Kondo impurities on the boundaries. Integrability of the models is established by relating the Hamiltonians 
to one parameter families of commuting transfer matrices. This is achieved through solving the reflection equations 
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for non-c- number solutions. Finally we solve the model by means of the algebraic Bethe ansatz method and derive 
the Bethe ansatz equations. 



II. GRADED REFLECTION EQUATION ALGEBRA AND TRANSFER MATRIX 

In this section, we give a brief review about the Z2-graded boundary quantum inverse scattering method. To 
begin, let V be a finite-dimensional linear superspace. Let R G End(V (g> V) be a solution to the Z 2 -graded quantum 
Yang-Baxter equation (QYBE) 

Ri 2 (ui - u 2 )Ri 3 (u 1 - u 3 )R 23 (u 2 - u 3 ) = R 23 (u 2 - u 3 )Ri 3 (ux - u 3 )R 12 (ux - u 2 ). (II. 1) 

Here Rjk{u) denotes the matrix on V ® V ® V acting on the j-th and fc-th superspaces and as an identity on the 
remaining superspace. The variables u±, u 2 and u 3 are spectral parameters. The tensor product should be understood 
in the graded sense, that is the multiplication rule for any homogeneous elements x, y, x' , y' 6 EndV is given by 

(x <g> y)(x' ® y') = (xx' g> yy'), (11.2) 

where [x] stands for the Z2-grading of the element x. Let P be the Z2-graded permutation operator in V ® V . Then 
P{x®y) = (— l)[*Ibf]y ®as, Vx,y G V and R 21 (u) = Pi 2 Ri 2 {u)P 12 . 
We form the monodromy matrix T(u) for a L-site lattice chain by 

T(u) = R 0L (u)---R m (u). (II.3) 

Indeed, one may show that T(u) generates a representation of the Z2-graded quantum Yang-Baxter algebra, 

Ri2(ui - u 2 ) T (ui) T (u 2 ) =T (u 2 ) T (ui)Ri 2 (u-l - u 2 ), (II.4) 
where for notational convenience we have 

T(u)=T 13 (u), T(u)=T 23 {u), 

and the subscript 3 now labels the quantum superspace W = V® L . 

In order to describe integrable Kondo impurities in strongly correlated electronic models with open boundary 
conditions, we need to introduce an appropriate Z 2 -graded reflection equation (RE) algebra. We introduce the 
associative superalgebras 71 and T + defined by the i?-matrix and the relations 

1 2 2 1 

Ri 2 (ui - u 2 ) T- {ui)R 2 i(ui + u 2 ) T- (u 2 ) =T- (u 2 )R 12 (ui + u 2 ) T- (ui)R 21 (ui - u 2 ), 

r£ tst2 (- Ul + u 2 ) rf 1 (ui)iei 2 (-u a - u 2 + v ) m 

=r; st2 (u 2 )r 21 (- Ui -u 2 + v ) rf 1 ( Ui )r$ lst *( Ul - u 2 ). (11.5) 

Here 77 is the so-called crossing parameter and sti stands for the supertransposition taken in the i-th space, whereas 
isti is the inverse operation of s tj. On e of the important steps towards correctly formulating the Z2-graded case is to 
introduce in the second RE in ( 1I.10| ) the inverse operation of the supertransposition. In our cases, the R- matrices 
enjoy the unitarity property, 

i?i 2 (ui - u 2 )R 2 i(-ui + u 2 ) = 1, (II. 6) 

and the crossing unitarity 

R s £ ( Ul - u^Rll 1 + u 2 + r?) = 1. (II.7) 
One can obtain a class of realizations of the superalgebras T + and 71 by choosing T± (u) to be the form 

71(u) = T.{u)f.{u)Tl\-ul Tf(u) = T?(u)f«(u) (T+\-u)) st , (II.8) 

with 
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T_(u) =R OM (u)---Roi(u), T + (u) = R 0L (u)---R 0<M +i(u), f ± (u) = K±(u), (11.9) 

where M is any index between 1 and L, and if±(u), called boundary if -matrices, are representations of 7±. In the 
following, without loss of generality, we shall choose M = L so that T + (u) = K + (u). 

The if -matrices if ± (it) satisfy the same relations as 7± (u) , respectively. That is the if -matrices obey the following 
REs 

1 2 2 1 

Rl2(ui - U 2 ) if- {ui)R 2 l(ui + U2) if- (112) =K- (u 2 )Ri2(ui + 112) if- (ui)i?2i(wi - u 2 ), 

R£ lst2 (~ Ul + uj,) iff 1 ( Ul )R 12 (- Ul -u 2 + r)) iff* 2 (ua) 
2 1 
=if; si2 {u 2 )R 2l (-u l ~u 2 +ri) if;* 1 {utjRfi lst2 (~u 1 +u 2 ), (11.10) 

Following Sklyanin's approach |30|, one defines the boundary transfer matrix t(u) as 

t(u) = str{K + {u)T_(u)) = str (K + {u)T{u)K_{u)T- l (-u)) . (11.11) 

Then it can be shown that (2l]] 

[t(ui),t(u2)]=0. (11.12) 

Although many attempts have been made to find c-number boundary if matrices, which may be referred to as 
the fundamental representation, it is no doubt very interesting to search for non-c-number if matrices arising as 
representations in some Hilbert spaces, which may be interpreted as impurity spaces. 



III. INTEGRABLE NON-C-NUMBER BOUNDARY if -MATRICES AND KONDO IMPURITIES IN THE 

ONE-DIMENSIONAL T — J MODEL 



Let Cj a and Cj. a denote creation and annihilation operators for conduction electrons with spin a at site j, satisfying 



the anti-commutation relations given by {ej CT , Cj. a >} 



6ij5 



where i, j = 1, 2, • • • , L and c, a' =T, 4- We consider 



the following type of Hamiltonians describing two magnetic impurities coupled to open t-J chains 



L-1 



L-1 



j=i,<? i=i 
J a Si • S Q + V a ni + J b S L ■ S b + V b n L . 



(IILl) 



Above, the projector V 



nLa- 



-n^riji) ensures that double electron occupancies of sites are forbidden; J a , V a (a 



a, b) are the Kondo coupling constants and the impurity scalar potentials respectively; Sj as usual is the vector spin 
operator for the conduction electrons at site j; S a (a — a, b) are the local moments with spin-s located at the left and 
right ends of the system respectively; nju is the number density operator n,j a — c^ a Cj a , rij — nji + riji- 
For the choices 



t=l, J = 2, V 



1 

2' 



(m.2) 



it has been shown in refs. ]22|-[24| that the bulk Hamiltonian acquires an underlying supersymmetry algebra given 
by <7^(2|1) in the minimal representation. Throughout we will refer to this case as the supersymmetric t-J model. 
Integrability of this model on a closed chain with periodic boundary conditions was established independently in works 
by by Essler and Korepin [^H and Foerster and Karowski [^6) by showing that the model can be constructed using 
the QISM. Furthermore, open chain integrability with appropriate boundary conditions was shown in refs. [p7[-p9[. 
It is quite interesting to note that although the introduction of integr ablc i mpurities we propose below spoils the 
supersymmetry, there still remains su(2) symmetry in the Hamiltonian (III.l) which mai ntains conservation of total 
spin and electron number. We will establish the quantum integrability of the Hamiltonian (III.l) for the special choice 
of the model parameters (ill. 2) and 



(c a - s)(c a + s + 1) 



V n = 



s(s + 1) 



(c a - s)(c a + s + 1) ' 



(III.3) 
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This is achieved by showing that this model can be derived from the QISM. Our result is consistent with the appli- 
cability of the coordinate Bethe ansatz method discussed in Eg) , 

Another choice of couplings which leads to an integrable model on the closed periodic chain is given by 



t = 1, 



J 



-2, V 



3 
2' 



(III.4) 



as shown by Schlottmann p2| corresponding to an su(3) invariant solution of the Yang-Baxter equation. In this case 
we can introduce integrable Kondo impurities on the boundary by choosing 



(c a + s)(c a - s - 1) 



l-s(s + l) 



(c a + s)(c a 



1) 



(III.5) 



Below, we describe how the aforementioned integrable cases are obtained. 

Let us recall that the local Hamiltonian of the supcrsymmetric t — J model is derived from an i?-matrix satisfying 



\ 






, (III.6) 
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where u is the spectral parameter and 



P = £(-l)W e < 



is the form of the Z2-graded permutation operator in accordance with the rule (II.2). We chose to adopt the Z2-grading 

[1] = [2] = 1, [3] =0on the indices labelling the basis vectors. 

We now solve ( H-ICQ for K + {u) and K-(u). For the quantum R- matrix ( III.6| ), One may check that the matrix 
K- (u) given by 



A-(u) B-(u) 
K_{u) = \ C-{u) 
1 



(III.7) 



where 



A-{u) 
B-(u) 
C-(u) 
D-(u) 



Cgjcg + l)-u 2 + u- s a (s a + 1) + 2uS z a 

(c a + U- Sa)(c a + U + S a + 1) 

2wS^ 

(c a + u- s a )(c a + u + s a + 1) ' 
2uS+ 



(c a + u- s a )(c a + u + s a + 1) 
c a (c a + 1) - u 2 + u - s a (s a + 1) 



2uS* 



(c a + u- s a )(c a +11 + s a + 1) 



(Illi 



satisfies the first equation of ( |II.10 )(for the details, see the Appendix). Here — S x ± iS y . The matrix K+(u) can 
be obtained from the isomorphism of the superalgebras 71 and T + . Indeed, given a solution K- of the first equation 
of (pT0|), then K+(u) defined by 



Kf(u) = K_{-u+-) 



(III.9) 
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i s a so lution of the second equation of (|II.10| ) . The proof follows from some alg ebra ic com putat ions upon substituting 
( IH.9 ) into the second equation of ( |II.1C| ) and making use of the properties flll.6| ) and ( II. 7 ) of the i?-matrix with 
rj = 1. Therefore, one may choose the boundary matrix K + (u) as 



A+(u) B+(u) 
K+(u) = | C+(u) D+(u) 
1 



(111.10) 



where 



A + (u) 
B+{u) 
C+(u) 
D+(u) 



s b (s b + l) + (2u-l)S| 



(Cb + U + Sb)(cb + 

(2« - l)S b " 



s b - 1) 



(c b + u + s b )(c b + u - Sb - 1) ' 

(2m-1)S+ 
(c b + u + s b )(c b + u - Sb - 1)' 
cg-u 2 ~s b (s b + l)-(2u-l)Sg 

(Cfc + M + Sb)(cb + U — Sb — 1) 



(III.ll) 



It can be shown that the Hamiltonian (III. 1 ) is related to the logarithmic derivative of the transfer matrix t(u) 
with respect to the spectral parameter u at u = (up to an additive chemical potential term) 



L-l ! 

H = J2 h iJ+i + o K '- (°) 



atr A' + (0)g £ o 



(111.12) 



with 



h = —PR(u). 
au 



For this case we obtain (III.l) subject to the constraints (III. 2, [II. 3). This implies that this model admits an infinite 
number of conserved currents thus assuring integrability. 

The second choice of integrable couplings results from use of an i?-matrix obtained by imposing Z2-grading to the 
fundamental su(3) i?-matrix and which reads 

\ 



R{u) 



u - 1 





























— u 





-1 























u 











1 











-1 





— u 





























— u — 
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u 





1 











-1 
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1 


































u + 



(111.13) 



V 

where again u is the sp ectral parameter and we adopt the same choice for the Z2-grading of the basis states as before. 

We now solve ( |IL10| ) for K + {u) and K_{u) for this R- matrix ( 111.13 ). One may check that the matrix K_(u) given 
by 



A-(u) B-{u) 
K_{u) = ( C-{u) D_(u) 
1 



(III. 14) 



where 



A-(u) 
B-(u) 



Ca(Cg - 1) - U 2 - U - S a (Sg + 1) - 2llS Z a 

(c a -u + s a )(c a - u - s a - 1) 

2wS^ ^ 

(c Q - u + s a )(c a - u - s a - 1) ' 
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CL(u) = - 
D-(u) 



(c a -It + s a )(c a -u - s a - 1) ' 
Cg(c a - 1) - u 2 - u - s a (s a + 1) + 2uSl 

(c a - U + S a )(c a - U- S a - 1) 



satisfies the first equation of (11.10). For this case K + (u) defined by 

K$(u)=K_(-u-\) 

is a solution of the second equation of ( ILlCj ), since the crossing parameter r] 
boundary matrix K + {u) as 

A+(u) B+(u) 
K+{u) = | C+(u) D + {u) 
1 



where 



A+(u) 
B+(u) 
C+(u) 
D + (u) 



2u - s b (s b + 1) 



(2w + 3)Sf 



(c b - u + s b - l)(c& 
(2 U + 3)S b - 



u - s b 



(c — u + Sb- l)(cb - u - s b - 2) : 



(2w + 3)S 6 



(cfc - u + s b — l)(cb — u — Sb - 2) 
eg - u 2 - 2u - s b {s b + 1) - 1 + (2u + 3)S^ 



(c b - u + s b - l)(c 6 



s 6 -2) 



(III. 15) 
(111.16) 

-3. Therefore, we may choose the 

(111.17) 



(111.18) 



As before, the Hamiltonian in this case is defi ned as the logarithmic deriv ative of the transfer matrix at zero spectral 
parameter giving the same general form (III. 1 ) subject to the constraints (III. 4, III.5| ). 

Before concluding this subsection,we would like to point out that the quantum integrability of other two cases 
corresponding to the same t but with negation of all J's and V's follows from the transformation — » (— l) J 'd, 

Cj (-iy Cj . 



IV. THE BETHE ANSATZ SOLUTIONS 



Having established the quantum integrability of the models, let us first dia gona li ze th e Hamiltonian ( III. 1 ) by means 
of the algebraic Bethe ansatz method |]30| , ^l| for the choice of couplings ( [II. 2 , III. 3 ). We introduce the 'doubled' 
monodromy matrix U (u) , 



A n (u) A 12 {u) Bi(u) 
U{u) =T(u)K-(u)T(u) = \ A 2 i(u) A 22 (u) B 2 (u) 

Ci{u) C 2 {u) V{u) 



(IV.l) 



where T(u) — T 1 (— u). Substituting into the graded reflection equation ( 11.10] ) , we may draw the following commu- 
tation relations, 



Ad(ui)C c (u 2 ) = 



(ill — U 2 — l)(u\ + U 2 ) . ,,.1, , \4hn t \ * r \ 

) -^r( Ul +u 2 + l)f h r( Ul ~ u 2 y c h d C e (u 2 )A gt ( Ul ) 

(Ul - U 2 )(Ul + u 2 + 1) J 

( U1 + ^ + i)tr + i) (aB + i) r(a - + 'Kfe(*w«J + 

5 _||_ y ,( 2 „ l + 1 )g Ca (« 1 )A e («=), 



2?(wi)C b (u 2 ) 



(ui -u 2 - l){ui + u 2 ) 



(ui - u 2 )(ui + u 2 + 1) 

— — — r C d (ui)Adb(u2) 

ui + u 2 + 1 



C b {u 2 )V{ Ul ) + 



2u 2 



(ui - u 2 )(2u 2 + 1) 



Cb(ui)V(u 2 ) 



(IV.2) 



(IV.3) 



G 



Here Abd(u) = Abd(u) + 2m +i &bdD{u) and the matrix r(u), which in turn satisfies the quantum Yang-Baxter equation, 
takes the form, 

rg(u) = l, r M = __L- ) rg(u) = ^,(6^d,6,d = l,2). (IV.4) 
Choosing the Bethe state as 

|n) = c ll K)---c lJV ( U7V )|*)^ 1 - lN , (IV.5) 

with being the pseudovacuum, and applying the transfer matrix t(u) to the state |f2), we have t(u)|0) = A(u)|f2), 
with the eigenvalue, 

AC 1 - 2 " Z 1 ( Cb + " ~ Sb ^ ( C b + U + S b + I) U + 1 L -|^|- (U + - Mj - 1) 

W ~~ 2u + 1 (c b + it + s 6 ) (c 6 + u - s b - 1) u - 1/ Al ( u _ Uj )(u + u . + i) 

2u , U 2 L T^r {u + — Uj — 1) 



_^ n ; u+uj i! u :" j :?i ^)( U; k», ( iv. 6) 



2ti + 1 u 2 — 1 j-J-^ (u — Uj)(u + uj + 1) 
provided the parameters {%} satisfy 

2Uj - 1 (q, + Mj - Sj,) (Cfc + Mj +£6 + 1) , Uj + l _ r U 

zmj (C& + Uj + Sb) (Cb + Uj - Sb — 1) Uj 



(IV.7) 



Here A^(u;{ui}) is the eigenvalue of the transfer matrix t^(u) for the reduced problem, which arises out of the 



r(u) matrices from the first term in the right hand side of ( |IV.2|) , with the reduced boundary K matrices K± (u) as, 

\C W {u) D^{u) )' 



where 



A^\u) = 


c 2 - 


u 2 - 


Sa{Sa + 1) + 


(2u- 


+■ i)S2 


(c a 


+ u 


- S a ){c a + U 


+ S a 


+ i) 


B ( l\u) = 






(2«+l)S- 






(c a 4 


- u — 


Sa)(Ca + U + 


Sa + 


i)' 


C W (u) = 






(2u + l)S+ 






(c a 4 


- u — 


Sa){Ca + U + 


Sa + 


i)' 


D W (u) = 


c 2 - 

'-"a 


u 2 - 


Sa(Sa + l) - 


{2u- 


f i)s; 


(c a 


+ u 


- S a )(c a + U 


+ S a 


+ i) 



and 



where 



(IV.9) 



4 l) w 


9 9 


- S b (Sb + 1) + 


(2u 


-i)Sf 


(c b + 


u + s b )(c b + u 


- s b 


-1) 






(2u - 1)S 6 " 






(c b + u 


+ s b )(c b + u — 


Sb - 


1)' 






(2u - 1)S+ 






(Cb + u 


+ s b )(c b +u — 


Sb - 


1)' 






- s b (s b + 1) - 


(2u 


-i)Sf 


(c b + 


u + s b )(c b + u 


- s b 


-1) 



(IV.ll) 
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Here K^\u), the boundary K matrices after the first nesting, follows from the relations, 



A dd (u)\^) 
A db (u)\$) 



2u 



2u + 1 
2u 



2u 



-)(- 



u 



1" u 2 -l 



2u+ l" aD y '~"~' " v "'" uv u 2 - 1 

Indeed, applying the monodromy matrix T(u) and its "adjoint" T{u) to the pseudovacuum, we have 

T«m(u)|*> = u L \q), T33(U)|*) = {u+ 1) L |*}, 
T 3d (u)|*}^0, T d6 (u)|*} = 0, T d3 {u)\*) = 0, 

f dd (u)\V) = (--JL_)*|*>, T33(tt)|*> = 



(IV.12) 



r 



r 3 ,i(ti)|*> # o, T d6 (u)|*) = o, f d3 ( u )|*) 



: 0. 



1' 



(IV.13) 



Then we have 



2?(u)|#) 
B d (u)\9) 



1 



(- 



M 2 - I 



)l 

i J fC_( U ) dh |*), 



v It- 1 

0, C d (u)|*>^0, 



f 



) L (K_(u) dd - 



1 



2u+ 1 



)l*> 



1 



2u + 1 



(- 



1 , 



1 



l*>- 



(IV. 14) 



In our calculation, use of the following relations has also been made 

(2u+l)T 13 (u)f 31 (u)+T 11 (u)f 11 +T 12 (u)f 21 (u) = -(2u - l)f 31 T 13 (u) + f 32 T 23 (u) + f 33 T 33 (u), 
(2u+ l)Ti 3 (u)f 32 (u) + Tn(u)f 12 + T 12 (u)f 22 {u) = -2uf 32 T 13 (u), 

T 2 i(u)f n + T 22 {u)f 21 {u) + (2u + l)f 23 T 31 (u) = ~2uf 31 {u)T 23 {u), 
T 21 (u)f 12 {u) + T 22 (u)f 22 + (2u + l)T 23 (u)f 31 (u) = f 3 iT 13 (u) - (2u - l)f 32 T 23 (u) + T 33 T 33 (u), 

which come from a variant of the (graded) Yang-Baxter algebra (1X4) with the R matrix ( III.(i| ), 



(IV. 15) 



T (u)R(2u) f (u) =f (u)R(2u) T (u). 



(IV.16) 



Implementing the change u — * u + i with respect to the original problem, one may check that these boundary K 
matrices satisfy the reflection equations for the reduced problem. After some algebra, the reduced transfer matrix 
may be recognized as that for the iV-site inhomogeneous XXX spin-^ open chain with two impurities of 
arbitrary spin on the boundaries, which may be diagonalized following Ref. p0|. Here we merely give the final result, 



A* 1 ) («;{«,•}) = 



(cb+u - s b ) (c b + u + s b + 1) 



(c b + u + s b ) (cb + U- s b 



n 



1) - LJ - u + c a 

a—a,o 



-{ 



2u-l 
2u 



M 

n 



(u - v m + f )(u + v m + |) 

j (U - V m + \)(U + V m - i) 



2u + 1 -p-r (u - c a + s a ) {u + c a + s a ) -pr (u — u 3 ){u + uj + 1) {u — v m — \)(u + v n 



2 ii 



a — a , b 



s a ) (u + c a - s a ) AJ- (u - Uj - i)(u + Uj ) JLi^ ( u _ Vm + + Vm _ i) 



n 



2 > 



l\ J ' 



(IV. 17) 



provided the parameters {v m } satisfy 



1\ N , 



n 



(v r , 



M 

n 

fc=i 



(Vrn -Vk~ l)(^m + Vk ~ 2) 
(«m - «fc + l)(Um + «fc) 



(IV.18) 
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After a shift of the parameters Uj — ► Uj — \,v m — > u m + ^,the Bethe ansatz equations (fV.7) and (IV. 18) may be 
rewritten as follows 

/%__2\2L TT M 3 + c a + 5 a + 2 _ TT U j ~ U ™ ' 2 U J ~"~ m ' 2 
I j 11 „._,. _ c i ll 



Uj 2 a=ab U j c a s " 2 m =i U J Vm 2 U J + u m 2 



TT Vrn-C a - S a V m + C a - S a yr (v m " Uj ~ |) (« m + Uj - |) _ yj (^m ~ ^fc ~ 1) Kn + ^fe ~ jj 

Q ia 6 Vm ~ Ca + Sa Vm + C " + S « jJi ( Wm ~~ U l + I) ( Wm + U + f ) fli ( Wm _ U fe + !) ( U ™ + «fe + 1) : 

with the corresponding energy eigenvalue E of the model 

N 



E 



E^r- ( IV - 19 ) 



1 U 7 



It should be pointed out that when s a = ^,the above results reduce to those obtained in Ref. [BO l , which in turn 
provides an algebraic interpretation for the applicability of the coordinate B ethe ansat z method |18|j. 

We now perform the algebraic Bethe ansatz procedure for the couplings ( III. 4 , III.5| ). We introduce the 'doubled' 
monodromy matrix U (u) 

/ A u (u) Ai 2 (u) Bi(u) \ 
U{u)=T{u)K_{u)f{u)= A 2 i{u) A 22 (u) B 2 {u) , (IV.20) 

\ d(«) C 2 (u) V{u) ) 

where T(u) = T~ l (—u). Substituting into the reflection equation ( [I.lOj) , we find the following commutation relations 
A b d{ui)C c (u 2 ) = — U2 + s)^ U ] + U2 , %p- r(ui + u 2 + l)f h r{ui - u 2 )f d C e {u 2 )A gi {ui) + 



2(«i + l) , ^gb 



(ui - u 2 )(2ui + 1) 



r(2«i + l)SC 9 («i)Ao(«a), (IV.21) 



V( Ul )C b (u 2 ) = j^- — U2 + U2 \ c b (u 2 )V( Ul ) + ^| — C h {ux)V(u 2 ) 

(Ul — U2)(Ul + U2 + 1) (Ui — u 2 )(2u 2 + 1) 

1 CaiuMdb^). (IV.22) 



Ml +712 + 1 



Here A b d{u) = A b d{u) + 2 u+i $bdV{ u ) and the matrix r(u), which in turn satisfies the quantum Yang-Baxter equation, 
takes the form, 

rg(u) = l, rg = -J-, r b d t(u) = ^—,(b^d,b,d= 1,2). (IV.23) 

U + 1 71+1 

Choosing the Bethe state |f2) as 

|fi) =C tl (u 1 )---C lN (u N )\V)F^-^, (IV.24) 

with 1^} being the pseudovacuum, and applying the transfer matrix t(u) to the state \Q), we have r(u)|fi) = A(u)|f2), 
with the eigenvalue 

A(u) — 2,1 + 3 (Cb-U + S b ) (Cfc - 7i - fife - 1) 71 + 1 L -pr (u + Uj)(u - 77j - 1) 

2u + 1 (cb — u + s;, — 1) (cb — u — Sfe — 2) u — 1 (u — Uj)(u + Uj + 1) 

2u U 2 L T^r (it + Uj + 2){U — Uj + 1) 



v 7i 2 - r 11 (u-u^fu + u^ + ll v ; 



2ii + 1 7i 2 — 1 (u — Uj)(u + Uj + 1) 

provided the parameters {u^} satisfy 



9 



M 



2uj + 3 (Cb-Uj + Sb) {Cb - Uj — 86 — 1) ^JLj + 1 ^2Z, TT (Uj + Ui)(uj — Ui — 1) 



2uj + 2 (ci, - Uj + Sb — 1) (ch - Uj - s& - 2) 



L _-^ («j + Ui + 2)(uj — u.i + 1) 



■A (1) («j;{«i})- ( IV - 26 ) 



Here A^(u;{ui}) is the eigenvalue of the transfer matrix t^(u) for the reduced problem, which arises out of the 
r(u) matrices from the first term in the right hand side of (IV. 21), with the reduced boundary K matrices K±\u) as 



K ( l\u) = 



A W (u) B W {u) 



C a) (u) D^'(u) 



(IV. 27) 



where 



A W {u) 



B W (u) 
C^\u) 
D { }\u) 



eg - u 2 - 2u - s a (s a + 1) - 1 - (2u + l)Sg 

(c a - U + S a )(c a ~ U — S a — 1) 

(2^+l)S- 

(c n - u + s a )(c a - U - S a - 1) ' 

(2^+l)S+ 

(c a - u + s Q )(c Q - u - s a - 1) ' 

c\ - u 2 - 2u - Sajsg + 1) - 1 + (2u + l)Sg 

(C a ~ U + S a ){c a - U — 8 a — 1) 



(IV.28) 



and 



4 i )(u) = MV ; M 4>) 



(IV.29) 



where 



4 x) ( U ) 



cf(«) 



eg - m 2 - 2m - s b (s 6 + 1) - 1 - (2u + 3)Sg 
(e& - U + Sb -l)(cb — U — Sb — 2) 

(2u + 3)S~ 

(cb - u + Sb - l)(cb - u - Sb - 2) ' 

(2tt + 3)S+ 

(c b - u + s b - l)(c b - u — Sb - 2) ' 
c 2 - u 2 - 2u - s 6 (s b + 1) - 1 + (2u + 3)S| 



[Cb — U + Sb — l)(Cfe — U — Sb - 2) 



(IV.30) 



Here K_ (u) ,the boundary K matrices after the first nesting, follows from the relations, 



Add{u)\$) 
Adb{u)YS>) 



1 



2u+l" u — 1 



1*) 



u 2 - 1 • 



(IV.31) 



Indeed, applying the monodromy matrix T(u) and its "adjoint" T(u) to the pseudovacuum,we have 

T dd (u)\V) = u L |*), T 33 (u)|tf) = (u + 1) L |*>, 
T 3 d(«)|*> + 0, T db (u)\V) = 0, T d3 (u)\V) = 0, 

f M (u)\*) = (--^-r) L \*), fssWi*) = (-4 ± r) i l*) I 

u z — 1 u z — 1 

r 3d (u)|*) #0, f db (u)\9)=0, f <a («)|*) = 1 



(IV.32) 



Then we have 
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2>(u)|#) = (- 



M+l, flT , 



u-V 

B d (u)\*) = 0, C d (u)\V)?0, 

,2 



Adb{u)\V) = (- 
Add(u)\V) = (- 



- 1 



'iT-(u)«»l*)> 



1 



-)l*> 



1 



1 ^ // - I 



u 2 - 1' y " w "" 2u+ 1" ' ' 2u 
In our calculation, use have also been made of the following relations, 

(2it + l)Ti 3 (u)T 3 i(u) + Tii(«)Th + Ti 2 (u)T 2 i(u) = (2it + l)T 31 Ti 3 (u) + f 32 T 23 {u) + f 33 T 33 (u), 
(2u + l)T 13 (u)f 32 (u) + Tu(u)f 12 + T 12 (u)f 22 (u) = 2uf 32 T 13 {u) 1 

T 2 i(u)f u + T 22 {u)f 21 {u) + {2u + l)f 23 T 31 {u) = 2uf 31 {u)T 23 {u) 1 
T 21 {u)f 12 (u) + T 22 (u)f 22 + (2u + l)T 23 (u)f 31 (u) = f 31 T 13 (u) + (2u + l)f 32 T 23 (u) + f 33 T 33 (u) 



which come from a variant of the (graded) Yang-Baxter algebra ( [1.4 ) with the quantum R matrix ( III. 13 ), 



(IV.33) 



(IV.34) 



i t t i 

T (u)R(2u) T (u) =T (u)R(2u) T (u). 



(IV.35) 



Implementing the change u — > u + \ with respect to the original problem, one may check that these boundary K 
matrices satisfy the reflection equations for the reduced problem. After some algebra, the reduced transfer matrix 
T^'{u) may be recognized as that for the iV-site inhomogeneous XXX spin-^ open chain with two impurities of 
arbitrary spin on the boundaries, which may be diagonalized following Rcf. pOjj. Here we merely give the final result, 



A (1) («;{«j}) = 



(c& - u + s b ) (c b — u-s b — l) 
(cb - u + Sb - 1) (cb —u — Sb — 2) 



n 



U + C a + Sr 



a—a : b 



M 



f 2u + 3 -^r (u-v m - + v m + 1) 
2m + 2 (u-v m )(u + v m + 2) 



(u - c a - s a + 1) (u + c a - s a + 1) 



2u+ 1 yr 

2u + 2 - 1 (u — c a + s a + 1) (u + 

a— a,o ' 



N 

n 



(u — tij)(u + Uj + 1) 
j= Jj (m-Uj + 1)(u+ u j +2) 

provided the parameters {v m } satisfy 

JV 



n 



(u - Um + l)(u + Vm + 3) 
(U - V m )(u + V m + 2) 



(«m + C Q - S a + l)(« m - C 



a =a 6 + Ca + Sa + ~ c a + s a + 1) ^ («m ~ «j + l)(«m + «7 + 2 ) 



1) tt («m -Uj)( v m + «j + 1) 



}• 



M 

n 

k=i 



(IV.36) 



(Vm ~ Vk ~ l)(^m + ^fc + 1) 
(«m - Wfe + l)(fm + Wfc + 3) ' 



After a shift of the parameters Uj 
rewritten as follows 



(IV.37) 

l,the Bethe ansatz equations ( IV.26| ) and (IV.37) may be 



i) 2i n 



a— a.o 17 



1 N 

m 



(iij - u. ( - l)(uj + Ui - 1) 



n 

a—a.b 



N 

n 



1) 



M 

n 



2 "3 



2 j 



Dm + 



+ 3) («m + u j + \) 



n 

fe=l 



(Vm ~Vk - I) (Vm + ^fc - 1) 
(«m - "fc + 1) («m + «fc + 1) : 



with the corresponding energy eigenvalue E of the model 



1 u ? 



(IV.38) 
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V. CONCLUSION 



In this paper, we have studied an integrable Kondo problem describing two impurities coupled to the one-dimensional 
t — J open chain for specific couplings. The quantum integrability of the system follows from the fact that the 
Hamiltonian may be derived from a one-parameter family of commuting transfer matrices. Moreover, the Bethe 
Ansatz equations are obtained by means of the algebraic Bethe ansatz approach. It should be emphasized that the 
boundary K matrices found here are highly nontrivial, since they can not be factorized into the product of a c-number 
K matrix and the local monodromy matrices. However, it is still possible to introduce a "singular" local monodromy 
matrix L{u) and express the boundary K matrix K_(u) as 



K-{u) = L(u)LT l {-u), (V.l) 



where 



— c — 1 — S z — S~ 

L(u)=\ -S+ Q u-c-l + S* ] , (V.2) 



which constitutes a realization of the (graded) Yang-Baxter algebra fln.4| ) when e tends to 0. The implication of such 
a singular factorization deserves further investigation. Indeed, this implies that integrable Kondo impurities discussed 
here appear to be, in some sense, related to a singular realization of the (graded) Yang-Baxter algebra, which in turn 
reflects a hidden six-vertex XXX symmetry in the original quantum i?-matrix. Therefore, one may expect that the 
formalism presented here may be applied to other physically interesting strongly correlated electron systems, such as 
the supersymmetric extended Hubbard model |52| and the supersymmetric U model Q. However, our construction 
is not applicable to the one-dimensional Hubbard model and the one-dimensional Bariev model, although a hidden 
six- vertex symmetry occurs in these two physically interesting strongly correlated electron systems ]3^| . Moreover, the 
singularity of the local monodromy matrix ( |V.2| ) implies that we can not apply it to construct a closed t — J chain 
interacting with integrable Kondo-like impurities. This is different from the conclusion by Zvyagin and Schlottmann 
fll| , j35f , who claimed that integrable magnetic impurities exist in the closed t — J and Hubbard chains. 

As shown in , one can put integrable Kondo impurities on the boundaries of the 5-function interaction electron 
gas. Obviously, there should be no problem in applying our construction to this model. Another question is to extend 
the present construction to the g-deformed case. This will lead us to integrable anisotropic Kondo impurities coupled 
to the q-deformed version of the t — J open chain. 

In concluding, we would like to point out that it will be interesting to carry out the calculation of thermodynamic 
equilibrium properties of the model under consideration, based on the Bethe ansatz equations presented here. Es- 
pecially, it is desirable to calculate the finite-size spectrum analytically, which, together with the predictions of the 
boundary conformal field thery, will allow us to draw various critical behaviour properties. Also, our construction 
may shed new light on a long-standing problem about the quantum integrability of the conventional Kondo model by 
QISM, given it has been solved using the coordinate Bethe ansatz for a long time MM- 
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APPENDIX A: DERIVATION OF THE NON-C-NUMBER BOUNDARY K-MATRICES 



In this appendix, we sketch the procedure of solving the (Z2-graded) RE for K-(u) . To describe integrable Kondo 
impurites coupled with the one-dimensional t — J open chain, it is reasonable to assume that 



K-{u) 



A{u) B{u) 
Ciu) D{u) 
1 



(Al) 



Throughout, we have omitted all the subscrips for brevity, refl ecting that the fermionic degrees of fr eedom do not 
occur, as it should be for a magnetic impurity. For the i?-matrix (III. 6), one may get from the RE (11.10) 33 functional 
equations, of which 11 are identities. After some algebraic analysis, together with the su(2) symmetry, we may assume 
that 
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A(u) = a(u) + [3(u)S z , B(u) = j3(u)S- , 
C(u)=p(u)S + , D(u) = a(u) -(3(u)S z . 

There are two equations automatically satisfied, leaving only 20 equations left to be solved 

A(ui)B(u 2 ) + B{ux)D{u 2 ) = A(u 2 )B( Ul ) + B(u 2 )D( Ul ), 
C{u x )A{u 2 ) + D( Ul )C{u 2 ) = C(u 2 )A( Ul ) + D(u 2 )C( Ul ), 



(A.2) 



u_(A(m)B(u 2 ) + B{m 
u_(A(u 2 )B( Ul ) + B(u 2 
u^{C{ui)A{u 2 ) + D(u x 
u-{C{u 2 )A{ui) + D(u 2 
u-(A(ui)A(w) +B(m 
u-(A{u 2 )A(ui) +B(u 2 
u-{C{ Ul )B{u 2 ) + D( Ul 
U-(C(u 2 )B( Ul ) + D(u 2 
u_((u+ - l)B(ui)D(u 2 
«-((«+ - l)C{ux)A{u 2 
- l)A{u 2 )B{u x 
u_((u + -l)D(u 2 )C(«i 



D(u 2 ))=u+(B( Ul )-B(u 2 )), 
D(u 1 ))=u+(B{u 1 )-B(u 2 )), 
C(u 2 )) = u + (C(u 1 )-C(u 2 )), 
C{ Ul ))=u+(C( Ul )-C(u 2 )), 
C(u 2 ) - I) = u+(A( Ul ) - A(u 2 )), 
C( Ul ) - 1) = u+(A(ui) - A(u 2 )), 
D{u 2 ) - 1) = u+(D(ui) - D(u 2 )), 
D(ui) - 1) = u+(D(ui) - D(u 2 )), 

- A{ux)B{u 2 )) = u+{{u- - l)D{u 2 )B{ Ul ) 

- D{u x )C{u 2 )) = u+((u_ - l)A(ua)C(tti) - 



- D(ui)B(u 2 )), 
A( Ul )C(u 2 )), 
■ B(u 2 )A(ui)), 



- B{u 2 )D{u x )) = u+((u_ - l)B{ Ul )A(u 2 ) 

- C{u 2 )A{ Ul )) = u+((«_ - l)C{ Ul )D{u 2 ) + C(u 2 )D( Ul )), 

u_(A( Ul )A(u 2 ) - (u+ - l)B( Ul )C(u 2 ) + (u + - l)C(« a )B(«i) - D(u 2 )D( Ul )) = u+(D(u 2 )A{ Ul ) - D{ Ul )A{u 2 )) 
u-(A(u 2 )A( Ul ) - (u+ - l)B(u 2 )C( Ul ) + (u+ - l)C( Ul )B(u 2 ) - D( Ul )D(u 2 )) = u + {A( Ul )D(u 2 ) - A(u 2 )D( Ul )). 
A( Ul )B(u 2 ) + u+u^D( Ul )B(u 2 ) + («_ - l)A(u 2 )B( Ul ) = (u+ - l){B{ Ul )D{u 2 ) + (u_ - l)B{u 2 )D{ Ul )), 
D( Ul )C(u 2 ) + u + u-A( Ul )C(u 2 ) + (u_ - l)D(u 2 )C( Ul ) = (u+ - l)(i7(«i)A(«a) + («- - l)C(tta)A(«i)), 
B(u 2 )£>(iii) + u+u_B(m 2 )A(ui) + (tt_ - l)fl(ni)£>(u 2 ) = (u+ - l)(A(u 2 )£(ui) + (u- - l)A(u t )B(u 2 )), 
C(u 2 )A(ui)+u+u-C(u 2 )D(ui) + («_ - l)C(«i)i4(u2) = («+ - l)(D(u2)C(«i) + («- - l)D( Ul )C(u 2 )), 

with u+ = Ui + u 2 , U- = U\ — u 2 . Substituting ( |A.2| ) into these equations, we find that all these equations are reduced 
to the following three equations 

u + {a{u\) - a(u 2 )) = u_(— 1 + a(ui)a(u 2 ) + s(s + l)/3(ui)(3(u 2 )), 
u + (P(ui) - (3(u 2 )) = u-(a(ui)0(u 2 ) + a(u 2 )0(Ui) - (3(ui)0(u 2 )), 
u+(a(u 2 )(3(ui) — a{ui)(3(u 2 )) — u-{a(ux)(3(u 2 ) + a(u 2 )f3(ui)) + u-(u + - l)(3(ui)0(u 2 )). 



(A.3) 



Taking the limit U\ 



u 2 , these equations become 

da(u) 

du 
d/3(u) 

du 
da(u) 



^(-l + a(u) 2 + S ( S + l)P(u) 2 ), 
2u 



a{u) 



d(3{u) 



■/3(u) 



1 

2u 
1 



(2a(u)/3(u) 



f3(uf), 



du du 
Solving the first two equations, we have 

(cic 2 -u 2 )(2s + l) + (c 2 



= — (2a(u)f3(u) + (2u - 
2u 



a(u) 



ci)u 



(2s + l)(ci-u)(c2-u) 



I3{u) 



l)(3(u) 2 ). 



2(c 2 - c\)u 



{2s + l)(ci -u)(c 2 -u) 1 



(A.4) 



(A.5) 



where c\ and c 2 are integration constants. Substituting these results into the third equation in (A.4), we may establish 
a relation between cf. c 2 = ci + 2s + 1. This is nothing but the non-c- number boundary K matrix ( III. 7 ) (after a 
redefinition of the constant: ci ^ — c — s — 1). 

A similar construction also works for the quantum R matrix ( [11.13 ). 
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